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CALCULATION OF SMALL ROCKET TRAJECTORIES 
USING ONLY FOCKET-TO - LAUK C HER RADIAL VELOCITY DATA 


by 

J. E. Jackson and S. M. Radicella* 

National Aeronautics and Space Administration 
Goddard Space Flight Center 
Greenbelt, Maryland 

A very accurate method for obtaining rocket trajectories 
is the radio Doppler tracking system, known as DOVAP (Ref. 1 
and 2). This system provides radial velocities (and by 
integration radial distances) between the rocket and the 
receiving site. In the DOVAP system, which was developed 
by the Ballistic Research Laboratory, a minimum of three 
vrell- separated ground stations are required. Three sets 
of radial distances are obtained and the rocket position 
is calculated by a triangulation technique. A simplified 
version of this system has been developed by J. C. Seddon 
(Ref. 3) in which only one station is required and in which 
a two-axis interferometer is used to supplement the radial 
velocity data. Since the Doppler measurement gives the 
radial distance, and the interferometer gives the direction 
cosines, their joint use provides the rocket position vector. 
It has been found, however, that a satisfactory interferometer 
requires a fairly elaborate instrumentation, and that its 
calibration is a very difficult procedure. Hence extremely 
skilled personnel are needed for its operation. 


*0n leave from the University of Tucuman, Argentina 



In the present paper a method is developed for calculating 
rocket trajectories using only the Doppler data at one receiving 
site. With this method it is possible to calculate small 
rocket traiectories above the drag region (Z-’TO km) with an 
accuracy of about 1 km. In practice only a few simple operations 
on a desk calculator are required in making use of this method. 

To justify the method, a detailed and lengthy discussion 
of its underlying theory is required. However the simplicity 
of the method is readily seen from the summary which is given 
after the theory. The calculations are based upon a measure - 
ment at the launching site of the rocket radial velocity R 

versus time using Doppler data at 73.6 Me. If desired a 

,t . 

counter can be used to determine R = \ R dt and this 

Jo 

information can also be recorded during the rocket flight. 

The method should be applicable to any small rocket 
traiectory with a peak altitude less than 200 kilometers; 
however since it has been checked mainly with typical Nike- 
Cajun data, the present report is based upon Nike-Cajun 
trajectories. Examination of Nike -Cajun performance data 
(Ref. 4) shows that, at the time of the second stage burnout, 
the horizontal range is of the order of 4 km. Assuming a 
maximum azimuth dispersion of the order of +30 degree, will 
place the "plane" of the trajectory at most 2 km away from 
the launch site. Since peak altitudes are typically in excess 
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of iGO kjt \ f i r.6 T^vdicil dist&ncts "R to tnt i ^ tii.c t 3 1 1 p fs mvpn 
with sufficient accuracy by: 

R 2 = X 2 + Z 2 

where : 


X - horizontal displacement 
Z = altitude 

Even if the Doppler sire is as such a? 4 km a wav from the 
plane of the trajectory, the approximation made in formula 
(1) introduces an error cf less than 0.1 a in the value of R 
at a radial distance of 100 km. For consi sfencv with "he 

system of coordinates introduced later in Section III, the 
Z and X directions are defined as the vertical direction and 
the direction of the velocity vector, respectively, at the 
peak of the trajectory. The origin of the coordinate system 
is the projection of the Doppler receiving sire upon the 
XZ plane. In the simplified analysis (Sections I and II) 
the earth curvature and the Coriolis effects are neglected. 

I . Calculation of Peak Time 

Differentiating equation (1) vita respect to time gives: 

RR = XX + ZZ (2) 
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For 


0 and 


the peak time (t ). 

o 


' Z = 


(R") = (XX) (3) 

z o t o 

For the accuracy desired, the Coriolis acceleration can be 

neglected and we can assume that X is a constant. The error 

introduced by this approximation in the final trajectory 

determination is shown in section III to be less than 0.5 km 

for a typical Kiks-Cajun flight. Applying equation (2) for 

two times (t - t) and (t + t) symmetrical with respect 
o o 

to the peak time t gives: 
r o ^ 


and. 


• 0 


(RR) . = (ZZ)„ + (X - Xt)X 


• • 


- (zz, t +t + (x t_ + xt)x 


and since 


(ZZ) 


t Q +i 


= - (ZZ) 


tft-t 


(RR) r r + (RR). , . = 2 X X 
L o L t o‘ ^ r o 


Making use of (3), this last equation can be written: 


(RR). „ + (RR). - 2(RR) 

l o l 


V 1 - 


or 


(RR) + (RR) = 2(RR) 

- 1 +t o 


(4) 


- 4 - 


■where the following simplifications in notation have been made 
in the subscripts: (t Q ) changed to (o) , (t - t) changed to 
(-t) and (t + t) changed to (+t). In this simplified 
notation X , Z . and R will represent values of X, Z and 

o’ o o 

R corresponding to the peak time. 

The calculation of peak time wbiich is based upon formula 

(4) will be explained with the aid of Figure 1. This figure 

shows a curve of RR as a function of time for the free fall 

* 

portion of a typical Kike -Cajun trajectory. The RR values 

were obtained from equation (2), using theoretical values 

of ZZ and XX for the trajectory defined by the parameters 

shown in Figure 1 and assuming no Coriolis effects. The 

* 

ordinates 4, B, and C correspond respectively to (RR) , 

(RR) and (RR) . Thus: 

— t *4“ t 

B + C = 2A 

It should be noted that in the illustration C is a negative 
quantity, as will always be the case for times (t Q + t) 
remote from peak. Let the ordinates A’, B' and C' correspond 
to times (t Q + A* t) , (t - t + a* t) and (t Q + t + A’t). 

This is equivalent to over-estimating the peak time. It is 
readily seen that: 

B’ + C'gi B + C 

but that A’ < A 
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Hence: 

B' + C' > 2A' , if peak tine estimate is too large. (5) 

Conversely for ordinates A", B" and C", symmetrical with 
respect to (t -A”t): 

B" + C” < 2A", if peak time estimate is too small. (5)' 


The ordinates A”, B" and C M are not shown on figure 1. 
A convenient form of equation (4) is as follows: 


(HR) + 
max 


(RR) = 2A 
mm 


= 2(RR) 

o 


(4)' 


The determination of the peak time t can therefore be 

o 

* 

made in the following manner. First, the function RR 

is calculated and plotted. The values of (RR) and 

max 

(RR) fflin are read from the graph and added. This gives a 

very good estimate of the quantity 2A and consequently of 

the time t , which is the abscissa corresponding to the 

ordinate A. One could also obtain this estimate directly 

from interpolation between tabulated values of RR as a function 

of time. The graph is therefore mostly an aid in visualizing 

the operation and in smoothing the experimental data points. 

The value of t^ obtained by this method is usually within 

a second of the correct value t . 

o 

The next step could be to select symmetrical times 
(t^ - t) and (t^ + t), such that corresponding ordinates B 
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- and C will fce clone to (kk) 


• PP \ 


* nen one c o:np s 


the v?]i« 


of (E 


C) 


with the value of 2k. where A is now the 


ordinate corresponding to t ' , end use the criterion developed 
above (Formulas (5) and (5)*) to determine whether the first 

estimate of t should be increased or decreased. Bv 

o 

observing the manner in which the difference (B + C) - 2A 

is altered in successive trials, one can determine very 

rapidly (by extrapolation or interpolation) the value of t 

with an accuracy of +0.25 second. It should be noted that 

this is a rapidly convergent process, and consequently that 

the initial estimate of t need not be particularly accurate. 

Thus it is possible to begin with the theoretical maximum 

possible value of t , namely the value t at which R is 

maximum (RR =0). In practice this estimate would be too 

great by 5 to 30 seconds. In general the method based upon 
• • 

(RR) and (RR) yields the answer wnth sufficient 
max ruin 

accuracy . 

It is also of interest to note that pr 


and RR 

max min 


occur at an altitude of no less than ^ (Z Q ),that is^ above 
the crag region for a normal bile -Cajun trajectory. The proof 
of this statement is as follows. The maximum and minimum 
values of RR occur when 


(RR) = Z 2 + ZZ + X 2 = 0 


/ 


or 


7 2 y 2 72 C2 

Z - — - ■+• — = — ( 1 + 0 -) where - = 2L_ 
-Z -Z -7 7.2 


— 1 7 

Using the approximation Z = g, Z = gt and h = 1 gt gives 


Z - h = 2h(l + a) 
o 


or 


h = 


3 + 2a 


shoving that the distance dropped h is less than , 

" 3 

II . Calculation of the Peak Altitude 

This calculation is based upon formulas (2) applied at 

a tine t + t : 
o 1 

(HR) = (XX) + (ZZ) 

1 1 1 


( 6 ) 


The left term of equation (6) is an experimentally determined 

quantity. 

(RR) 

From equation (3): X = £ (7) 


and from (6): 


(RR) = X(X + Xt ) + (Z - h )Z 
1 o 1 oil 

Where h^ is the distance fallen from peak during the time 
interval t^. A method for calculating and h^ is given 
in the appendix. From (7) and (8): 


( 8 ) 
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0 


(RR) - (RR) 
1 < 


(RR) 


tj - (Z o - h )Zj 


or: 


(RR). - (RR) 

i < 


(R ^n 


R 2 - Z 


- t - (Z - h )Z = 0 
2 1 o 1 


(9) 


Equation (9) can be re-written as a third degree equation: 

A. Z 3 + B, Z 2 + C Z + D = 0 
i O 1 o 1 r» 


(10) 


Where : 


h = z , 


p = 

n 


(RR) - (RR) 
o 1 

- R 2 A 
o 1 


h Z 

1 1 


D = I (RR) - (RR) 1 P/ - t (RR) Z + R 2 h Z 
1 u 1 o J o o oil 


Then, Z^ is the solution of the equation: 

f(z) = A 1 z 3 + B 1 z 2 + C z + D = 0 (11) 

The function f(z), for the sane trajectory as was used in 
deriving figure 1, and for t = - 105 seconds is shown in 
figure 2, with coefficients based upon the correct peak time 
(£t - 0) and also based upon r 1 second errors in peak time 
(At = +1 , and /It = -1). It was found necessary to use four 
significant figures on R and five significant figures on R 
(Calculated points shown as circles) to define the f(z) 
function. This accuracy on R represents the very best that 
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can be expected from a Doppler measurement, however, in 
practice R can at best be determined to three significant 
figures. The f(z) function was therefore re-calculated for 
the At = 0 curve using the best accuracy which can be achieved 
experimentally and the calculated values were shown as crosses. 
It is seen that in practice it should be quite difficult to 
define the f(z) function, using experimental data. It is 
also seen that for At = -1 sec there are two solutions 

(Z = 168.4 km and Z = 173.9 km'), and that for At = +1 sec 

o o 

there ane no solutions. Even for At = 0 the cubic (11) has 
two positive solutions very close to each other, one of which 
is the correct answer. A criterion for selecting the correct 
solution, can be derived from the sign of the tangent to the 
curve representing equation (11). It can be shown that the 
sign of the slope is, at the true solution, the same as the 
sign of the time t ^ or, in other words, it is negative for 
ascent values and positive for descent values. The above 
discussion suggests the undesirability of using equation (11) 
directly for obtaining Z^ unless extremely accurate values of 
RR are available. 

It is of interest to note that, when the correct peak 

time is used, the minimum of the function f(z) occurs for a 

value of z close to Z , and that this minimum can be deter- 

o 

mined more accurately (by solving f’(z) = 0) than the f(z) 


10 



ru.iccicn. fos^i-vati on suggests that a fairl** good vaii 
or Z can be obtained by solving f • (z) = 0. A study of the 
derivative f omul a showed that 


i_ i i \J. L «J 


tv t. ■ — _i_ I i 


by proceeding as shown in the following analysis. 

A very accurate second degree equation (whose solution 
requires loss precise experimental values) can be obtained 
by assuming that 


X = XT 

o 


wnere i is the tree fall time corresponding to the peak 
^ti'ttice Z i.e. T is defined by the following equation: 

Z = \ g (1 + - T 2 )T 2 
o 2 o 

in which the term ~ T 2 is a correction for the variation of 
with altitude. The deviation of equation (13) is given in 
the appendix. The quantity X t which appeared in equation 
can then be evaluated as follows: 

From (12) 


i 2 = 


o X o g o° w2 ) 


X 


T7 


2Z 


-.2 

* 9 

V- =1 o 


-z, 


2Z g (1 + 2 - tf)t 
o fc o 1 


g (1 + - T ) 


and 


•2 Z, X" 

X t = - -i C 

2Z 


1 + 


1 + 2 a t J 


( 12 ) 

(13) 

g 

( 8 ) 


11 


( 14 ) 


* 2 

X t = 


z, u x") 

fi- 

zz 


where 


X = L ± . r Jr . =l + a. (T 2 - 2 t 2 ) . 

1 + 2 a t 2 

For the example used in the text and for t =-105 seconds, 

\ = 1.002 

Equation (8) then becomes; 


(RR) = XX + X 2 t + (Z - h )Z 
1 o loll 


= (RR) - Jl 


Z, X i 


2Z 


-2 + (Z - h. )Z. 

o 11 


Changing (15) to a quadratic in Z^ gives: 


2 Z 2 Z - 2 

[(RR) 

- (RR) + h 

Z ] 

|z o - z 

X X 2 

o 1 

1 

o 

1 J 

° 1 

o 

2 Z 2 + 7 Z 2 - 2 

" (RR) ! - 

! 

(RR) o . J 

z 

- 7 z 2 

- * X 


o 2+^ 


2 f (RR), - (RR) 1 

__ [ - Q + ‘ J Z 


R 2 = 0 


o 2+7 ° 


The solution of equation (16) is: 


Z = p’ + 

o 1 V 2A 


+ p; 2 


(15) 


(16) 


(17) 
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where 



(RR), - (RR) h, 

i Q + 1 

(2 + ^)Z 1 2+* 


Since * is extremely close to unity equation (17) can be 
written: 



(18) 


where 


(RE) , - (RR) h 

i o + 1 

3Z, 3 


If is unity, equation (16) can be written 


3 Z 2 Z - 2 f(RR) - (RR) + h Z 1 - R 2 Z = 0 
° 1 <- 1 o 1 o i 


and using the notation of equation (10) 


3 A Z 2 + 2B Z + C = 0 
1 o 1 o 1 

This is the derivative of equation (10). 

Thus the cubic of equation (10) has a minimum occurring 

for a value of Z verv close to the correct answer. In 

o 

other words the cubic has two roots very close in value, 
as seen in Fig. 2. The approximation made in equation (12) 
introduces a small error in Equation (15). Actually 

• *7 

X Q = XT + £ -which is equivalent to. changing in equation (14) 
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.to (X - £)% cr (X" - 2 C X + 
o o c 


. ihus equation (14) 


be c ome s ; 


• 9 -Z. * 

X Z t = — 1 


2Z 


A I rj “ ^ ^ j 

-* h ' W + -3 J 


x rt X- 
o o 


* ? 
-v/ -Z, X z 


1. Q X f 1 -2 £ ] 

2Z L X 


f rJ 2C 

This introduces a small correction (= ^-) on the value of 

o 

-■ in ecuation (17). It can be shown that the value of Z 

o 

obtained by using * = 1 (instead of K - 1 - _) is too large 

X 

by an amount given by: ° 

,2C Xt 
e = (_) (-£-) 

X ° 


? ' 

For a Nike -Cajun trajectory^ — , is always less than 1 per cent. 

X 

For the trajectory used °to derive Figure 1, a one per cent 
change in * introduced a change in Z q of 0.1 km. Thus formula 
(18) is a very good approximation. 

O m * 

In formula (18) R*-, (RR) , and (RR) are experimentally 

o’ 1 o 

determined quantities; h^ and Z^ are calculated using the 

A(t) and B(t) functions given in the appendix for the time 

t^ corresponding to (RR)^ and for the value of g corresponding 

to Z . Since Z is unknown, it is necessary to calculate 
° o 

(18) using first an estimated value of g which will give an 
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estimate of Z . This estimate of Z can be used to re- 


o o 

calculate g and a better determination of Z is then possible. 

o 

In practice only tvo attempts are sufficient since a good 
first estimate of g is given by: 


g(R ) < g(Z ) < g( ypv 2 - (RR) t ) 

o o v o op 

in which t is the peak time. This criterion for estimating 

P 

Z is based upon the fact that Z is less than F but greater 

O . O O 

than yj - (RK)^t . The second part of the preceding statement 
follows from the following relations: 

Z 2 = R 2 - X 2 ; X X = (RR) and X > . 

o o O O O O t 

p 

III . Errors Due to Approximations bade in the Theory 
Discussion cf errors 

For more rigorous treatment considering the effects 
of earth curvature and rotation, it is convenient to use 
coordinate axis X', Y ^ and Z’ defined as follows. The Z* 
axis is the vertical direction at peak and the X' axis is 
in the direction of the horizontal velocity V at peak. 

The Y’ axis is chosen such that X', Y' > Z' form a right-handed 
system of coordinates. The time t is measured from peak 
and it is positive for increasing X'. From this definition, 
the rocket coordinates at peak satisfy rhe following relations: 
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X' = 0 


Y' = 0 


Z’ = Z 
o o 


X* = V 

o H 


Y’ = 0 
o 


Z' = 0 
o 


Due to the earth’s curvature, the Doppler station will 

+x2 

be at a distance — 2_ below the X'Y’ olane. where X is 

2R e ‘ ° 

defined as previously and R g is the radius of the earth. 

The relationship between this system of coordinates and the 

one used previously is seen in Fig. (3) to be as follows; 


X = X + X' 
o 


Y = d + Y' 


Z = Z' + 


2R 


Thus the main effect of the earth's curvature is to introduce 

X 2 

the term in the expression for Z. If X is less than 

2R e x 2 ° 

50 km, the term — Q— is less than 0.25 km and usually 

2R x 2 

negligible; however if X is 100 km, the term .. P. is about 

° ^e 

0.8 km and it may be desirable to take it into consideration. 

The effect of the earth’s rotation is to introduce a 
Coriolis acceleration which is defined in terms of the earth's 
rotation vector x.. This vector has a magnitude of 
0.7272 x 10"^ radians/sec and it is directed towards the 
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I 


North Star. At a given location it has a component ■= x sir 
in the vertical direction and a ^ = x cos 0 in the North 

iv 

direction, where 0 is the latitude (0 is positive in the 
Nor therm hemisphere and negative in the Southern hemisphere). 
The Coriolis acceleration a is: 


• G 


a = - 2 x X R 


i x X 
X 


j x Y 
J Y 


k a Z 

Z 


and its components are: 


a = 2 (x Y - x Z) 
X Z Y 


2 (x Z - x X) 
X Z 

2 (x X - x y) 
Y X 


The effect of the Coriolis acceleration can be illustrated 

by using an eastward trajectory, in which case x = x and 

Y N 

x = 0. Furthermore by assuming that the trajectory is at 
X 

middle latitude in the Northern hemisphere, x^ and x are 
positive quantities, having comparable magnitudes. This 
case will also represent a typical trajectory at Wallops 
Island. The Coriolis velocities are then: 


V = + 2\ (x_ Y - x Z)dt 

X jo Z N 
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\ m * a ■“ ^ ^ 

' = x A at 

V I 7 


40 


v ^ 2t * . X dt 

V * ^ \ 

i o *- 


These expressions c o -n he further simplified considering that 

• • • 

Z >;> Y and X constant giving: 


r = _2* (Z - Z ) = + 2 * h 
X NO N 


V = - 2 x X t « - 2 x X' 
Y Z Z 


V„ = + 2 x X t = 2 * X' 
Z N N 


The corrections on X, Y and Z due to the Ccric-1: 




A X 


- f 

n; d 


g t 2 dt 


^ 2 
3 


AX= ± 


Ay = -2 a l Xt dt 


£ 


AY = - A z Xt 4 


AY = - * X't 


AZ = 2 X { Xt dt 


Az = x x’t 

N 


For a typical small rocket trajectory the region of interest 
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corresponds to values of t between -150 seconds and +150 
seconds, representing the altitudes between the peak and 100 km 
below the peak. Taking the maximum values of t = + 150 seconds. 


h = 100 km, X’ = 50 km, and jl 


Jl v 


= 0.5 x 10"^ gives: 


z 41 

AX = +0.5 km for +150 seconds and -0.5 km for -150 seconds. 

Ay = -0.375 km for t = +150 andAZ = 0.375 km for t = +150. 

It is seen that even the maximum values of AX.AY^andAZ 

are quite small. However, if desired, a correction can be 

made for the Coriolis effect, if the approximate direction 

of the plane of the trajectory is known. An accuracy of 

+10 degrees for this information is adequate to estimate the 

X’ , Y’ , Z* components of * and hence calculate the Coriolis 

terms. From the above discussion it is seen that AX, AY, 

AZ, V , V , V are all expressed in terms of time and flight 
x y z 

parameters derived from the approximate analysis fox a plane 
non-rotating earth. The approximate analysis could therefore 
be refined by letting: 


X = X + X' t + AX 
o o 


* * 

X = X + V 
o X 


X^ 

Z = Z’ + 2R~ + A2: Z = I g ol B(t ) + V z 

e 

and introducing corresponding correction terms in Equations 
4 and 18. Thus Equation (4) would become: 
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( 19 ) 


(RR) + (RR) - 4X * h = 2(RR) 

+t -t o N 

For the trajectory used in Fig. 1, the term 4X a ..h is equal to 

oh 

o 

0.8 km /sec which increased the peak time value by 0.28 second. 

Similarly p in Equation (18) should be changed to p* , where: 

1 1 



( 20 ) 


where A Z and V are the Coriolis terms. To evaluate the 
Z 

correction introduced by changing p to p', the data in Figure 1 

1 1 

were re-calculated for an East trajectory and with the 

corresponding Coriolis terms. Starting with these more exact 

values of RR and using equation (18), it was found that the 

calculated value of Z was too low by 0.4 km. However using 

o 

the value of p* given in formula (20) the error in Z was 
1 o 

less than 0.1 km. In making this calculation the exact peak 
time was used, to insure that the results would show only the 
effect of the Coriolis acceleration. 

In summary, the Coriolis terms introduce for the example 
given a correction of 0.28 second in peak time and 0.4 km 
in peak altitude. Thus a slight improvement in accuracy can 
be achieved if desired. However even with this correction 
it is believed that there will be an uncertainty of about 
0.5 km in the final trajectory determination in view of the 
approximations made in formula (1), (18), and in the A(t) 
and B(t) functions. 
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- 7 Summary of Method 

1. From the Doppler data calculate and plot a graph 
of RR similar to that shorn in Figure 1. 

2. The peak time t is obtained by calculating 

(RR) = % (RR) + (RR) ] 0 

o *» max min 

• • 

and reading the time at which RR = (RR) q . One should attempt 
to read this time to within an accuracy of 0.25 second. 

3. The peak altitude is obtained from: 


where ; 


p i ■ T + 


= p ■+ hP- + p 2 
1 V 3 P 1 

h, (pi) - (RR) 


(18) 


* 

In Equation (18) the quantity (RR)^ corresponds to a convenient 

time t of the order of -100 seconds, giving values of RR 

close to the maximum of the RR function, and h and Z are 

1 1 

given by ; 


h = g A(t ) 
1 Z 1 
o 


Z l‘ 8 Z * ( V’ 

o 

where the A(t) and B(t) functions arc 


at ed in the apoendix 


A preliminary calculation of (18) will be required using an 


estimated value of g based upon: 

Z o 


g (R ) " g (Z ) < g ( / R 2 - (RR) t ) 

O o Vo Op 
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. The horizontal distar-''^ ^ < r.pak and t’-e horizontal 


velocity are thu determined respectively from 


Y 
* » 

o 






and 

X 

o 

5. The complete trajectory can then be calculated 
versus time f r om : 

2= Z - g\Z ) A(t ) 
o o 

and 

X = X + X t 
o o 

6. For a slightly better accuracy a correction for 
Coriolis effects can be made using the method given in 
Section III. However this correction is not required if an 
accuracy of the order of 1 km is satisfactory for the values 
of X and Z versus time. 
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Appendix 1 


If Coriolis effects are negligible; free fall calculations 
for sinall rocket trajectories can be performed with a very 
good accuracy by means of the formulas; 

h = g ^ [l + a tn = g A(t) 
o 2 ■ ° 


h “ g D t 


1+2 



g 0 8(t) 


where: h = distance fallen from peak 

h = vertical velocity 

g Q = value of g at the peak altitude Z . 

t = time from peak 

7 2 

a = 2.42 x 10 -/ /sec for a typical small rocket 
trajectory. 

If desired the value of h and h(=z) can subsequently be 
corrected for Coriolis effect as indicated in section III 
of the report. For convenience the A(t) and E(t) functions 
have been tabulated for a representative value of c equal to 
2.42 x 10 /sec . The derivation of the formula' for the A(t) 
and B(t) functions is as follows; 


Let g 


(Z ) 
o 


(R + Z ) 
e o 


where g 

e 


g at surface of the earth 
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and R = earth’s radius. 


Then 


g(Z - h) = g 
o 


e (R o + Z - h) 
e o 


tJ 


(R 


+ Z R + Z J 


e o 


e o 


= g(Z Q ) + — — 


= 8<V + 


R^ + Z 
e o 


?g(V h 

R + Z 
e o 


) 


Thus 


g(Z Q - h) = g(Z Q ) + kh 


The value of k is relatively constant for values of Z q 
ranging from 100 to 200 km as can be seen from the following 

table: 


z 

o 

100 km 

150 km 

200 

km 

k 

2.94 x 10' 6 

2.87 x 10‘ 6 

X 

CD 

CM 

10" 6 


Thus the calculation 
integrals 


of h and h consists in evaluating the 


■£[■ 


h = \ I g(Z ) + kh| dt 
o 


'] 


24 


and 


f ; 

1q 


h dt 


Since the contribution of the kh term is small, an estimated 

n 

value of h given by h = %g,„ -.i is adequate for calculating 
the correction term kh. Thus g can be written : 


= g(Z ) [ 1 + %kt Z J 


The integration yields readily: 


B (Z o 5t [ 1+ 5 t2 ) 


h = %g(Z c )t Z [l + i t/j 

Thus the constant - in the expression for A(t) and B(t) is 


given by: 


12 


and the values of ~ for the 100 to 200 km altitude range vary 
between 2.45 x 10 ^ and 2.34 x 10 ' . A value of - of 2.42 x 10 ^ 
which corresponds to an altitude of 125 km is considered 
representative of the free-fall portion of most '^ike-Cajun 
trajectories and it was used m calculating tne aa w one i . - ) 
functions . 

It should be noted that if T is the total falling time 
from peak (assuming no drag), equation (21) gives: 
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f 


Z ~ "5 g i, t 

O o 

the calculation of Z . 

o 


which is the Equation (13) used in 
A plot of Equation (13) is given in 


Figure 4. 

Values of g(Z ) for Wallops Island are given on the 
o 

following table: 


E 

/I N 


0 

.0097986 

60 

.0096173 

70 

.0095874 

60 

.0095577 

90 

.0095281 

100 

.0094986 

110 

.0094693 

120 

.0094401 

130 

.0094110 

140 

.0093821 

150 

.0093533 

160 

.0093246 

170 

.0092961 

180 

.0092677 

190 

.0092394 

200 

.0092113 
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Table of A(t) and B(t) Functions 


t A(t ) h B(t) Z 


0 

0 

0 


5 

12.50 

5.00 


10 

50.00 

10.00 


15 

112.51 

15.00 


20 

200.02 

20.00 


25 

312.55 

25.01 


30 

450.10 

30.01 


35 

612.68 

35.02 


40 

800.31 

40.03 


45 

1013.00 

45.04 


50 

1250.76 

50.06 


55 

1513.60 

55.08 


60 

1801.57 

60.10 


65 

2114.65 

65.13 


70 

2452.92 

70.17 


75 

2816.33 

75.20 


80 

3204.96 

80.25 


85 

3618. 82 

85.30 


90 

4057 . 94 

90.35 


95 

4522.34 

95.42 


100 

5012.10 

100.48 


105 

5527.22 

105.56 


110 

6067.73 

110.64 
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** 

u 

A(t) 

n 

bit . ) 


115 

6633.66 


115.74 


120 

7225.13 


120.84 


125 

7842.03 


125.95 


130 

8484.56 


131.06 


135 

9152.69 


136.19 


140 

9846.45 


141.33 


145 

10566.01 


146.48 


150 

11311.31 


151.63 


155 

12082.29 


156.80 


160 

12879.36 


161.98 


165 

13702.21 


167.17 


170 

14551.01 


172.38 


175 

15425.97 


177.59 


180 

16327.02 


182.82 


185 

17254.19 


188.06 


190 

18207.76 


193.32 


195 

19187.42 


198.59 


200 

20193.60 


203.87 


205 

21226.20 


209.17 


210 

22285.27 


214.48 


‘ 215 

23371.13 


219.81 


220 

24483.38 


225.15 
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Figure 2. 


Figure 3. 
Figure 4. 
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The RP. function for a typical Nike-Cajun 
trajectory. 

Graphical solution of the cubic 
3 o 

f(z) = Az +Bz +Cz+D =0, using 
^ 111 

coefficients corresponding to a typical Kike -Cajun 
trajectory. The effect of a + 1 second error 
in peak time is also illustrated. 

Coordinate system used for trajectory calculations. 
Peak altitude Z vs total falling time T. 
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Based upon: 

X = .42 km/sec 


T 



TIME (seconds) 












= PEAK ALTITUDE (km) 







where 


(RR), - (HR) h, 

- 1 _Q + 1 


(2 + X)z 


1 


2+* 


Since ^ is extremely close to unity equation (17) can be 
written: 


z o = p i V-r + p i 


(18) 


where 


(RR), - (RR) h, 

p = 1 o + 1 


1 


3Z 


1 


If * is unity, equation (16) can be written 


2 * 
3 Z Z 
o 


- 2 |(RR) - (RR) + h Zj-R 2 Z =0 

1 L 1 o 1 1 J o 1 


and using the notation of equation (10) 


3 A Z 4 + 2 B Z + C = 0 
1 o 1 o 1 

This is the derivative of equation (10). 

Thus the cubic of equation (10) has a minimum occurring 
for a value of Z Q very close to the correct answer. In 
other words the cubic has two roots very close in value, 
as seen in Fig. 2. The approximation made in equation (12) 
introduces a small error in Equation (15). Actually 
X = XT + C which is equivalent to changing X^ in equation (14) 



